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Abstract 

The work developed in the paper concerns the multivariate fractional Brownian motion (m£Bm) 
viewed through the lens of the wavelet transform. After recalling some basic properties on the mfBm, 
we calculate the correlation structure of its wavelet transform. We particularly study the asymptotic 
behavior of the correlation, showing that if the analyzing wavelet has a sufficient number of null 
first order moments, the decomposition eliminates any possible long-range (inter)dependence. The 
cross-spectral density is also considered in a second part. Its existence is proved and its evaluation 
is performed using a von Bahr-Essen like representation of the function sign(t)jt|". The behavior of 
the cross-spectral density of the wavelet field at the zero frequency is also developed and confirms 
the results provided by the asymptotic analysis of the correlation. 
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1 Motivations and overlook 



The fractional Brownian (fBm) motion developed by [Mandelbrot and Nesd (jl968l ) has been extensively 
studied as the archetypal model of fractal signals. Many extensions have been proposed, trying to keep the 
simplicity of its definition while modeling more complex phenomena. For example, time-dependent Hurst 
exponent or rf-dimension al extensions have been introd uced which have respectively led to the multifrac 



tiona l Brownian motion ([Peltier and Levv-Veheli 119951) and the fractional Brownian sheet (JAvache et al. 



20021 ). Another extension consists in defining multivariate fractal processes. This extensio n is needed 
by niany applications ranging from economy to physics, pas sing by biology and neuroscience (JGil-Alana , 



2003 : 



Arianos and Carbone . 



2009 : 



Achard et al. 



2008 



20061 ). In all these disciplines, many modern sens- 



ing approaches allow to measure instantaneously different variables from complex phenomena. The need 
of multivariate signals models is crucial in order to model and understand these phenomena^ 



A multivariate extension of the fBm has been proposed recently in a very general setting bv lDidier and Pipiras 



(|201l[ ) with the help of operator self-similarity. The operator fractional Brownian motion is an operator 
self-similar Gaussian process with stationary increments. When the operator is diagonal, it is called 
the multivariate fract ional Brownian motion (mf Bm). We have particu larly studied this diagonal case 



(Amblard et al 



20111 ) , elaborating on the work of 



Lavancier et al. 



(|2009l ). In these works, the correlation 



structure of the mfBm has been studied. The increments process has also been studied, and we have 
exhibited its correlation and spectral structures, showing the possible existence of long-range dependence 
in correlation as well as in cross-correlation between components. 

In this paper, we pursue the study by wavelet analyzing the multivariate process. It is now well- 
accepted that wavelet analysis i s the r ight f r amework to dea l with monovar i ate fr a ctal signal s with 



stationary increme nts (JFlandrin 



Bardet et al. 



1988 



1992 



Tewfik and Kim 



1992 : 



Wornell 



1990; 



Fay et al. 



2009; 



20001 ). Wavelet transforms or decomposition provide a regularized differentiation of the 



processes, have a filter bank structure in perfect adequacy with 1// type of spectral behavior, may 
eliminate long-range dependence properties if the analyzing wavelet is properly chosen. Thus, studying 
multivariate fractal signals through the lens of the wavelet zoom is indeed natural, and we expect that it 
will be useful as well in revealing the interaction structure between the components of the mfBm. 

We thus concentrate on the correlation and on the spectral structure of the wavelet transform of 
the mfBm. The principal result of the paper is the explicit form of the second order statistics of this 
multivariate Gaussian random field. We study the asymptotic behavior of the cross-correlation function of 
the wavelet, and this allows us to prove that choosing a wavelet with at least two null first order moments 
eliminates any possible long dependence in the correlations. We pro ve the existence and ca lculate the 
cross-spectral density of the wavelet field. We thus extend the result of lKato and Masrvl (|l999l ) providing 
the existence of the spectral density of the wavelet tr ansform of the fBm. Th e proof uses a generalization 
of the von Bahr-Essen representation of |t|" (used bv lKato and Masrvl (|l999[ l) to the function sign(t)|t|". 
We also provide the behavior of the density at the zero frequency, corroborating the asymptotic result 
obtained in the time domain. We stress on the fact that our asymptotic results do not impose that the 
wavelet is a real function, nor that it has a compact support. 

The paper is organized as follows. In order to have a self-contained exposition, we recall some basics 
definition and results on the mfBm in Section 1. In Section 2, we set the wavelet analysis and look at 
some self-similarity properties inherited from the process. The full correlation structure is developed 
in Section 3 whereas the spectral counterpart is in Section 4. The needed generalized von Bahr-Essen 
representation is proved in the last section. 



2 Some facts on the multivariate fractional Brownian motion 



The p dimensional multivariate fractional Brownian motion (mfBm) x{t) is defined as the only Gaus- 
sian process having stationary increments and having components jointly self-similar with parameters 
{Hi, . . . ,Hp) S (0,1)^. The self-similarity property can be stated as follows: for any real A > 0, 
x{Xt) = X^ x{t) where H = dia.g{Hi, ... ,Hp) and A^ is intended in the matrix sense. The notation 
= stands for equality of all the finite-dimensional probability distributions. 

The cross-covariance structure induced by the multivariate self -similarity property and the stationarity 
of the increments has been firs t studied by 
to the Gaussian assumption, 
more simple way as follows. 



Lavancier et al. 



Amblard et al 



(|2009[ ). Theorem 2.1, without having recourse 



(|201l[ ) have parameterized this covariance structure in a 



Proposition 1 (Proposition 3, (|Amblard et al.Ll201ll )^ Let j,k e {l,...,p}, j ^ k, then there 
exists (7j > 0, (pjkifljk) G [— 1, 1] X M satisfying pjk = Pkj o,nd ijjk = ^Vkj, such that 



rjk{s,t) ■.= E[xj{s)xkit)] 



{Wjk{-S) + Wjk{t) - Wjk{t - s)} , 



(1) 



where the function Wjk{h) is defined by 



Wjk{h) 



{Pjk-Vjksign{h))\h\"'+"^ 
P]k\h\ +rijkh\og\h\ 



ifH, 
zfH, 



Hk = 1. 



(2) 



This result is also valid in the case j ~ k when setting pjj = 1 and rjjj = in ([2]); we thus recover 
the covariance structure of a monovariatc fBni. The parameter a^ is the variance of a fBm at time 1, 
Var(xj(l)), whereas pjk represents the instantaneous correlation between components j and k at time 1, 
i.e. E[xj{l)xk{i)]- The antisymmetric parameter r]jk is related to the time- reversibility property of the 
multivariate process. Indeed, t he mfBm is tirne reve rsible, i.e. x(t) ~ x{—t) in distribution for every t, if 



and only if rjjk = for all j, k (JAmblard et al 



20111) 



To ensure that the matrix given by ([T]) is the cross-covariance matrix of a process, the constraints 



Amblard et al 



i mpos ed on pjk and rjjk are not sufficient. A necessary and sufiicient condition, proved by 

(J201ll ). corresponds to the positive-definitencss of the Hermitian matrix with entries T{Hj + Hk + 1) x ^j^ 

where $,jk is defined by 



?. 



[jk 



Pjk Sin ( I (if^- + Hfe )) - ir^^k cos ( f {H, + Hk)) iiH,+Hk^l 
Pjk - i^Vjk if Hj + Hk^l, 



(3) 



where i = \/— T. For example, when p ~ 2 and i]jk = no condition is required for the correlation 
pi2 but when Hi = H2, and when Hi = 0.1 and H2 = 0.2 the correlation pi2 cannot exceed 0.514, 



see (Amblard et al 



2011) for more discussion. The problem of simul ation of such a process has been 



Amblard et al 



investigated in 

([1]) presents some examples in order to illustrate the process. 



(|201l[ ) using the algorithm developed by IChan and WoodI (jl999l ). Figure 




"^^'*>?^/X/,>JC^^^^ 




•V^^A-Av,^/"^ 




Figure 1: Examples of discretized sample paths of a time reversible {rjjk = 0) mfBm of length n — 1024, with 
p — 20 components. The Hurst exponents are equally spaced in [0.3,0.4] (left plot), [0.6,0.7] (middle plot) and 
[0.4, 0.8] (right plot). The correlation parameters are set to 0.7 (left and middle plot) and to 0.3 (right plot). The 
components are shifted artificially in the left plot for the sake of visibility. 



The covariance structure of the increment process (at lag 1) can be easily deduced from ([T]). When 
j — k, we obviously recover the covariance of the fractional Gaussian noise and the classical property 
that this process is short-memory when Hj < 1/2 and long-memory when Hj > 1 /2. In the multivariate 



Amblard et al 



( 2011 ) 



case, long-range (interdependence) may also appear in the cross-covariance. Indeed 
proved that for all j ^ k, the cross-covariance behaves asymptotically as |/i|^J+^'=^^ (up to a constant) 
meaning that the long- memory property arises as soon as Hj +Hk > 1 which can appear in three different 
situations: Hj — Hk = 1/2, Hj < 1/2 and H^ > I — Hj or Hj > 1/2 and Hk > 1/2. In those cases, 
some troubles may appear when it comes to infer parameters of the models from data. Indeed, long- 
range d^£endBncejnByJead_to_ve2;_sl^^ in many works. 



e.g. (JFlandrin 



1992 : 



Veitch and Abrv 



1999 : 



Bardet et al. 



2000 : 



CoeurioUvl 120011 ). having recourse to 



wavelet types of transformation is an elegant way to overcome the problem. Indeed, using wavelet types 
of transformation with a correctly chosen filter allows to extract the stationary part from the fBni and 
allows to "whiten" the increments. We describe such an approach in the following section. 

3 Wavelet Analysis: definition, stationarity and self-similarity 

The use of wavelet analysis in the understanding of the monovariate fractional Brownian motion, an d 
more generally for the study of f ractal signals, goes back to the early works of iFlandrinI ([1988, 



Tewfik and Kim 



19921) 



(|1992| ). IWornelll (|l990f ) to cite some but a few. 
The aim is now to analyze the multivariate fractional Brownian motion through the lens of the wavelet 
transform. We use the continuous wavelet transform here, but a similar analysis could be performed in the 
multiresolution framework using orthonormal wavelet bases. We will consider complex valued wavelets, 
not necessarily in the Hardy class, not necessarily with compact support. The hypothesis we impose on 
the wavelets will be detailed when needed. 

3.1 Definition and stationarity 

Let ip he a. complex wavelet function, let a > and fe G K and consider i^abi-) — a'~-^/^il;{{. — b)/a). Let 



di 



■ipa,b) , = / Xj{t)i>a,b(t)dt 



(4) 



the wavelet transform of the jth component of a multivariate fractional Brownian motion. ■0 denotes the 
complex conjugate of ip. In this section, we assume that conditions [CI] and [C2(2)] are satisfied, where: 



[CI] Admissibility condition: tp{t) G L^ and [V'('^)|^/|w| G L^, where ^ is the Fourier transform of ip. 



[C2(K)] f"i/)(t) e L^ for m = 0, 1, . . . , K. 



Condition [CI] ensures that V'(O) ~ and that j^%p{t)dt ~ 0. We note, as iKato and Masrvl (Il999l ). 
that under condition [C2(l)], the integral ([4]) is well-defined as a sample path integral and is a second-order 
random variable. This follows, since under [C2(l)] we have j^ \s\^\i'a.b{s)\ds < +oo,\/H E (0, 1). 

The aim of this section is to focus on the correlation between the wavelet transforms (at different 
scales and different times) of two components j and k of the multivariate fractional Brownian motion. 
The wavelet transform is a random field. It is clearly zero mean and Gaussian. We have for ai,a2 > 
and b,h eR 



^K, . 64-/1 ^L.b] = / rjk{tl,t2)^paub+h{tl)i^a2At2)dtidt2. 



ai,b+h a2,b 

Under [CI], and from ^ the last expression reduces to 

a.ak 



EK^,b+hd-a2,b] = TT- / Wjk{t2 - ti)iia^,b+h{tl)i'a2.b{i2)dtidt2 



Let r^(z;) := J^'ipai^i,j^h{u)il}a2.b{u + v)du be the correlation function between the two wavelets ipai,b+h 
and 'ipa2.b- Then we have 



cr,(Tfc 



EK.b+A2,b] --^ w,u{v)TAv)dv. (5) 

Note that [C2(2)] implies that for all the values of Hj and Hk, /jg \wjk{v)\ \T.^{v)\ dv < +oo. With two 
changes of variables, this may also be rewritten as 



^idiub+hda2,b\ = ^Tr- y/a^x Wjk{a2t2-aiti- h)tl;{ti)ip{t2)dtidt2. (6) 



If we interpret for fixed parameters oi and a2 , the quantity -B[(i;^ b+h'^a b\ ^^ ^^'^ cross-correlation between 
two signals, we observe that it depends only on the difference between the times at which it is evaluated 
(i.e. h). With the fact that the wavelet transform is a zero mean and Gaussian field, we conclude that 
d^^ and d^^ are jointly stationary signals. 

3.2 Self-similarity type property of the cross-wavelet transform 

The variance of the wavelet transform at similar scales for the fractional Brownian motion with Hurst 



parameter H exhibits some self-similarity. Indeed, it is proved in iFlandrinI (|l988l ) for example that for 
all 6 

Varidi,)^a^"+^x ( -^ [ \t2 ~ ti\^"^{h)ij{t2)dtidt, 

V ^ JR2 

We note here that the same behavior holds for the cross-wavelet variance. 



Proposition 1 Under the assumptions [CI] and fC2(2)], let h ^ and fix ai = 02 = a > 0. Then, 

EKA] = a^^+«^+^ (-^ z,.) and Corr[dl„di,] = -^, (7) 

where Zju ■= Jj^2 Wjk{t2 ~ ti)i'{ti)ilj{t2)dtidt2. 



Proof. Consider Equation (j6]). The result is obvious when Hj + Hk 7^ 1 since for any a > 0, 
Wjk{av) = a^^'^^''Wjk{v). Now, when Hj + Hk = 1, the result comes from Condition [CI] ensuring that 
/r2 Vjk X a{t2 - ti)log{a)tp{ti)jp{t2)dtidt2 = 0. ■ 

Let us observe that the instantaneous cross-wavelet correlation is independent of the scale. 

4 Cross-correlation structure of the wavelet transform of the 
mfBm 

For fixed scales, 01,02, we now specify the behavior of the cross- wavelet covariance (or correlation) as 
\h\ — >■ -f-oo. In particular, our aim is to exhibit the influence of the number of vanishing moments of the 
wavelet function on the asymptotic cross-wavelet covariance. Such a result needs the following assump- 
tion: 

[C3] The wavelet function has Af > 1 vanishing moments that is 

t"'ij{t)dt = for TO = 0, . . . , Af - 1 and / t'^^i;{t)dt ^ 0. 



We may now derive our result obtained as \h\ — > +cx). Let us first recall Landau notation: for two 
functions f{h) and g{h) defined on M, we denote by f{h)^g{h) as \h\ — > +00 (resp. f{h)^o{g{h)) and 
f{h)^0{g{h))) ii\\m\h\^+oof{h)/g{h) = 1 (resp. lim|/i|^+oo ./(/j)/.9(^) = and f{h)/g{h) is bounded 
for all h). 

Theorem 2 Assume [CI], [C2(2M+1)] and [C3] hold, then as \h\ -^ +00, we have 

where K{i;,M) := (^j,f )(aia2)^^ \J t'^''iP{t)dt\^ and 

,, . {P3k+rioksign{h)){"'+"'') ifHj+Hk^l 



rijkycSign(h) 
2M(2A/-1) 



ifH,+Hk = l. 



We notice that the equivalence stated has a meaning as soon as pj^ -\- rjjk sign(/i) 7^ when Hj + i/j, 7^ 1 
and as soon as -qjk 7^ when Hj + Hk — I. In the opposite cases, a careful look at the proof shows that 



the equivalence can be replaced by an upper-bound or more precisely E['i^ t+h'^a b\ ~ ''(1^1^ ^^^) 



Similarly to the fractional Brownian motion, Theorem [5] asserts that the higher M, the less corre- 
lated the wavelet transforms of the components j and k of the multivariate fractional Brownian motion. 
This has many implications. In particular, this suggests that estimating the instantaneous cross-wavelet 
correlation at a scale a may be efficiently done by using the empirical correlation since at scale a, dP^ ^^ , ^ 
and d'l ^ are not too much correlated if M is large. 



Proof. The proof is split into two cases. Before this, we denote by Dh ■— {(^i, ^2) G K^ : [02^2 — ai^il < 

\h\ 



„ } and we note in particular that 

V(ii,t2) eDh 



a2t2-^aiti I < 1 < 1 and sign(a2t2 - ai^i - M = -sign(/i). 



Case 1 . a:^ Hj+Hk ^ 1. 



02^2—0.1*1 



< 2. 



We assume here that pjk + rijksign{h) 7^ 0. Let us write E[d-'^ b+h'^a b\ ~ — -^^•y/aia2 x T with 



T 



(Pjk - Vjk sigii(o2t2 - aiti - h)\a2t2 - aih - h\°'ip{ti)'ip{t2)dtidt2 = Ti + T2 



and where Ti (resp. T2) corresponds to the integral on Dh (resp. M^ \ Dh). Let us first prove that 
|^|2M-ay^ ^ as \h\ -^ +00. Denoting c^ = \pjk\ + \Vjk\, we have (since 2M - a > 0) 

h 



\h\^^^-"\T2\ < c^ 



\a2t2 - aiti\°'\h 



2M- 



'■\Du 



1 



a2t2 — aiti 



\ip{tl)\\lP(t2)\dtidt2 



< 22*^-°3"c^ 



'■\D,, 



ia2t2-aitiY^'\^{ti)\\i;{t2)\dtidt2. 



The result is then obtained by using assumption [C2(2M)] and the dominated convergence theorem. Now, 
within the domain Dh, one may use the series expansion of (1 + x)" (for |a;| < 1). 

/* / f / \ *^ 

Ti = |/ir / (pjk-Vjk sign{a2t2- aiti- h)) i I ^-^—- — ^1 tp{ti)ili{t2)dtidt2 



\hnpjk+vik 



sign{h)) f (j2{-r 



1,^ a\ f a2t2 - aih 
£ [ h 



Ipitl)ij{t2)dtidt2 



where (") denotes the binomial coefficient {a){a — 1) . . . {a — i + 1)/^!. Decompose Ti into three terms 
(denoted by T{,T2 and Tg) corresponding to the 2M first terms of the series, the (2M + l)th term 



{i — 2M) and the remainder terms. Then, 

2M-1 



i=o ^ ^ -^ 



(a2<2-aiii) tp{ti)i){t2)dtidt2 



Under Assumption [C3], ip has M vanishing moments and therefore the previous expression reduces to 

2A/-1 



2~i / \ r 

f—0 \ / J^^\Dh 



(02^2 - aiii) Tp{ti)%l^{t2)dtidt2- 



Now, 



2A/-1 






■,2A/-f 



x2Af 



2\Dh 



(02*2 - aiii)^" |V'(ii)IIV'(i2)|rftirfi2 



Assumption [C2(2M)] and the dominated convergence theorem may be combined to prove that \h\'^^^ "T[ 
as |/i| — > +00. The term Tj is defined as 



As previously we obtain 






Pjk + rijk sign(/i) \2M/ J^2 



(02*2 - fliii) tp{ti)ip{t2)dtidt2 



e'iP{t)dt 



2M 



K{i;,M). 



Since T = Ti + r2 = T{ + T2 + T^ + T2, the proof will be completed if we manage to prove that 



_ P]k + Vjk sign(/i) 



t' e>2M+i 



a\ / 02*2 — flii] 



'4'{tl)lp{t2)dtidt2 



f ia2t2-a,t,r'^'(j2i-lY+' 
■^^1^ \e>o 



Ct \ t 02*2 — aiti 

2M + 1 



'Ip{tl)'lp(t2)dtidt2 



The binomial coefficient appearing in the last equation satisfies, with i' — t + 2M + 1 

|a(a-l)---(a-f + 1)1 



2(2 ^ a) • • • (f -I -a) 
< — — -^j since a < 2 

2(f-l)! 2 2 
£'! i' - i 



Recall that in Dh wc have \a2t2 — aiti|/|/i| < 1/2. The series in the previous integral then satisfies 



E 

l>0 



i + 2M+lJ V 



f a2t2 - aiti 



< 



< 



2M 
2 



il-S 



£>1 



^ + 2i\/ + 1 



02^2 — aiti 



2M 



r + Si^- 



2M + 1 



+ 21og(2)=:a 



M- 



Thus we obtain 



|^|2Af-a|y,| < %^ / |a2t2-aiti|''^+V(ii)IIV'(i2)Miidi2. 



Since by Assumption [C2(2M+1)], t^'^'+H'it) e i\ we have |/ip*^""|T;^| = 0(|/i|-i), whence the result 



Case 2 . H^ + H^ = 1. 

We assume here that rjjk 7^ 0. We take the same notation as previously. We first note that, under 
[CI], the term T can be rewritten as 



T 



Pjk\a2t2 - aiti - h\ +rijk{a2t2 - aih - h)log 



02^2 — aiti 



h 



ij{ti)llj{t2)dtidt2. 



We decompose T in Ti + T2 (as done in case 1). The proof that \h\'^'^'^'^T2 ^ as \h\ -^ 
similar arguments as in the case 1 and is therefore omitted. Now, the term Ti can be rewritten as 



follows 



'Vokh 



Dh 



1 - --2^1-^^\ log f 1 - ^^H!1_^) Ipit,m2)dt,dt2. 



Denote by Ti and T2 these two terms. Assumption [CI] leads to 

02^2 — aiti\ -T 



Ti 



-Pjk\h\ I ( 1 



-P]k 



/ \a2t2-aiti\ 



Tp{ti)^it2)dtidt2 

dtidt2. 



Then, we assert that 



\hr-'\T,\ < 2''^-'3\p,k 



02*2 — aiti 



, 2 A/ 



'■\Dh 



i^a2t2 - aihj \ip{h)U{t2)\dtidt2 ^ 



as \h\ — > +00. For the term T2, we may use the series expansion of log(l + x) (for |a;| < 1). We omit the 
details and leave the reader to verify that as \h\ — > +00 

2A/-1 

T2 ^ rjjkh I I 1 - 



B2 



02*2 - Oiti 



1 f a2t2 — aiti 



2M-1 



10 



-1 /'a2t2__aih 
2M V h 



2M\ 






-h 



1-2M 



Vjk 



f2M\ 



2M{2M - 1) V M ) 



(0102) 



M 



r'ijj{t)dt 



Hence, T^ |/.|i-2M ^ (_Z|^2gSW) ^(^^^), 

In this proof, Fubini's theorem and interchanges of integrals and (in)finite sums are widely used. All of 
these are justified by the absolute convergence of the different scries related to the expansions of (1 + x)" 
or log(l + x) for \x\ < 1 and Assumption [C2(2M+1)]. ■ 



FlandrinI (1988, 



1992). A rigo r ous p roof of the existence of this spectral density in the L^ sense 



5 Cross-spectral density of the wavelet transform of the mfBm 

In the case o f the fBm, the expression of the spectral density of the wavelet transform was provided 

by 

was obtained by iKato and Masrvl (J1999I ). On the basis of this work, our ambition is to provide the 
cross-spectral density between wavelet transforms (at different scales) of components j and k of the 
multivariate fractional Brownian motion. The idea is to obtain the following spectral representation for 
the cross-correlation 



EK.b+kdt,,] 



1 

2^ 



Si';,,,iuj)e'-'du;. 



Theorem 3 Under Assumptions [CI], [C2(M)] and [C3J (with M > 2), we derive the following asser- 
tions, 
(i) The cross-spectral density of the wavelet transforms of two components j and k exists and is given by 



Ca,H = ^MH^a,afer(H, + Ffe + i)OfeH 



V'(aia;)V'(a2w) 



(9) 



whe 



Gfc(w) 



p,k sin ( f (if, + Hfe )) + i 77,fe cos {^(Hj+Hk)) sign{uj) if H, + Hk ^ I 
Pjk + i^Vjk sign{uj) if Hj + ifj. = 1. 



(ii) We have for both cases, as a; — > 
iiii) Moreover, the coherence function between the two components j and k satisfies: 



2M-l-a 



CtaM 



\SitaM\ 



S'^,aA^)S^ataM 



lOfeMl' 



T{2Hj + l)T{2Hk + 1) 



V'(aia;)V'(a2w) 
■0(aia;)^(a2w) 



(10) 



11 



Before writing down the proof, let us give some comments. 

1. Item (ii) in Theorem [3] is the speetral analogue of Theorem [21 Indeed the behavior ot the cross- 
correlation at infinite lags is linked to the behavior of its Fourier transform at the zero frequency. 
We recover the fact that as soon as M > Hj + Hk + 1/2, the long-range interdependence is 
destroyed. The divergence of \oj\^^^^^'^^'' is compensated by the rapid decrease to zero of the 
Fourier transform of the wavelet. 



2. The interpretation of the coherence (fTO| is difficult here. Indeed, it is complex valued, a property 
which is not natural for a coherence. This comes from the fact that the quantities S^ ^^ (uj) are not 
power spectral densities but cross-spectral densities (cross-spectral density between two different 
scales of the wavelet transform of one signal). Thus, to interpret correctly the coherence, we should 
look at one scale only, in which case we recover the coherence evaluated in the usual spectral domain. 
And this result is logical since the usual coherence is independent of the frequency. 

3. Setting ai = a2 and j = A: in the expression of the cross-spectral density, we recover the usual 
result of the power spectral density at one scale of the wav elet transform of a scalar fBm. The proof 
proposed here is a natural extension of the proof found in iKato and Masrvl (jl999[ ) . 



4. The derivation of the analytic form of the cross-s pectral density is easy if we use generalized func- 
tions (or Schwartz distributions). Indeed, from iGel'fand and Shilovl (|l964[ ). we know that the 
Fourier transforms of |i;|" and |f|"sign(?j) are respectively given by — 2r(Q; -I- l)|a;|~"~^ sin(7ra/2) 
and — 2ir(a + l)|ci;|~"~^sign(a;) sin(7rQ;/2). Then, in the proof below, the calculation of T = 
^^Wjk{v)TTp{v)dv can be done using Parseval equality. However, the theoretical background re- 
quir ed and hidden in the calcul ation is far more involved than the basics we have used in the proof 
Csee berfand and Shilovl \m^). 



Proof, {i) We recall that under [CI] and [C2(2)], Equation © holds, that is E[dl^,^^d^^,,J = -^T 
with T :~ Jj^Wjk{v)Tti,{v)dv. Furthermore, note that the Fourier transforms of xjja^b and r^(u) exist and 
are equal respectively to y/atlj{auj)e^^'^^ and to 



Now, let us split the proof into two cases. 



(11) 



Case 1. a := Hj + Hu ^ I. 



When 7 = fc, at t h is ste p, iKato and Masrvl (|l999f ) have used the representation of |z;|" obtained by 
von Bahr and EsseenI (|l965[ ). We have obtained a similar representation for the function sign(w)|w|" for 
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a E (0,2) \ {1} (see Equations P^ and (ITil) in Lemma S]). We have by Fubini's theorem and under 
Assumption [C3] (with il/ > 2). 



^ = {pjk - Vjk s^gn{v))\v\°'T.^{v)dv 

r(a + i) /•,,,,_„_! 



■^1 " / [ PjkSm{T:a/2){l — cos{luv)) ~ rijkC0s{T:a/2)a[gn{uj){sm{ujv) — ga{ujv)) \ X 
' |w| I -pjk sm(7ra/2) I I - rjjk cos(7ra/2)sign(a;) I I I duj 



- 7, 


r(a + i) 


TT 


r(a + i) 



|w| " ^ I pjfc sin(7ra/2) + i jy^-fc cos(7ra/2)sign(cj) jg(aj)(iaj 
kr"~^Cifc('^)'7(w)(ia;. 

Note that the condition M > 2 is required for a > 1. For a < 1, A/ > 1 is a sufficient condition. These 
conditions allow us to show that the contributions jT^,(v)dv and Jjg^ga{uJv)T^(v)dv are equal to zero. 
Now, using ([TT|) we obtain 



^"■JR"' V ' 

By using Bochner's theorem, the proof will be done, if one proves that the function P(-) is integrable. 
Let us prove this last assertion. Under [C2(M)], t^il>{t) E L^ for fc = 0, . . . , M. Therefore, -0 is a M times 
continuous and difFercntiable function. Using a Taylor expansion 

A/-1 

^(w) = ^ Lj'=^('^')(w)+cj^^^(*^)(w) =w*^?A(*^)(w), with we [OAw,OVw], 

fc=0 

under [C2(M)]. And since '0^*^^ is continuous at zero, 'tp{uj) ^ oj^^f/;'^*^) (0) as w — > 0. Then as w — >■ 0: 

P(w) ~ 0,(w)|a^P''-'-"(aia2)^^|^(^^)(0)p. (12) 

As a consequence, for M > 2, P is continuous at zero and lim(^_j.o± P(a;) — 0. Therefore for e > 0, P is 
integrable on the interval [— £, e] as a continuous function on this interval. Finally (with c^ := \pjk \ + |?7fcj | ) , 

|w|>e ^ J\u}\>ais PI ^ ^ J|w|>a2£ I'^l ^ 

under [CI]. Hence, P{-) G L^ and Bochner's theorem may be applied. 
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Case 2. H, + Hi. = 1. 

We start with the representation of v log \v\ given by ([T5)) . 

w'ifel") = P'jkV'\+-qjkv\og\v\^ Mm pjk\vr + 'njkv\og\v\ 

a— >1 

__ J_ /■ 2p.jk{l - cos{ix}v)) - 7r?7j7,sign(u;) sin(a;t;) 

Now, we derive the computation of the term T := j^Wjk{vyr^{~v)dv, similarly as the previous case. 
Using dominated convergence theorem and Fubini's theorem, 

T = ± f( lun f 2P.-.(1 - cos(c..)) - 7r^,.sign(..) sin(c..) ^^\ 

= ^ hm f ( f 2ft-^(l - cos(c.z;)) - 7r^,.sign(.;) sin(c..) \ ^ 



2^o!^-i [-^''' [ 2 j - -'?.'=«ign(-) (^ ^j ) ) l-l— dc. 

1 



lim / |w| " (2pjfc +i7r7]jfcsign(a;))7/;(aia;)'0(a2w)e'" dw. 



From (fT2|) . j^l " ^■!/'(aiaj)-0(a2'^) is an integrable function for all a G (0,2). Therefore, the integral and 
the limit may be interchanged. Therefore, we obtain 

1 f Pjk+i^Vjksign{uj)Y, 77. . ,u;h 



and Bochner's theorem can be applied. 
(m) is derived from (|T2|) . ■ 



6 Bahr and Essen type representations for the functions sign(i;)|v|'^, 

■u" and v^ 

In 1965, von Bahr and Essen have obtained the following representation theorem for |w|" for a E (0, 2): 

^ r(. + l)sin(W2) r i^f^rf^. (13) 

The following lemma provides a similar representationfor sign(w)|w|", u" = w"lK+(f) andw" = (— i;)"1r- (w). 
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Lemma 4 Let a G (0, 2) \ {1} and let ga : ffi — > M the function which equals zero when a G (0, 1) and 
which is the identity function when a G (1,2), then we have 

■ / M la r(Q; + 1) cos(7ra/2) [ signjiv) (sin(^t)) - ggjujv)) 
s^gn{v)\v\ = J^ ^^^ do;, (14) 

„ T{a + 1) f sin (tt^ ) (1 - cos(ww)) + cos (tt^ ) sign{uj) (sin(u;t;) - ga{ujv)) 

T{a + 1) /* sin (tt^) (1 — cos(ojf )) — cos (tt^) sign{uj) (sin(oji') — g{ujv)) 



The representations of w" and v" are obtained from ([T51) and ([T?|) noticing that 
< = ^ (li^r + sign(w)|w|") and «_ == ^ (!«!" " sign(f)|w|") . 

Proof. Let a G (0, 1), then from ([T5|) and properties of the function F 

1 I la+i r(a + l) /■l-cos(ww) 
-\v\ = cos(7ra/2) / — . . , „ — dio. 

Since Jj^ |a;|^"^^| sin(ajv)| < +oo for a G (0,1), we can differentiate this integral with respect to v to 

obtain 

• / M |Q r(a + l) /■ sign(a;)sin(wz;) 
sign(ti) u = ^ ^cos(7ra/2) / — dio. (15) 

When a G (1, 2), then from ([T5| and properties of the function F 

, ,„ 1 F(a + 1) , , ,^.. f 1 — cos(ww) , 
a\v\"-^ ^ ^ ^ (-cos(7ra/2)) / —^ — '-dw. 



Since L |w| " ^| sin(wu) — ujv\duj < +oo for a G (1, 2), we can take the primitive of the last equation to 
get 

• ^ M |Q F(a + 1) /■ sin(a;u)/a; - w 
sign(i;)|w|" = cos(7ra/2) / — ^—-^ duj 

^ Jr \^\ 

^(" + ^).„c^ /r,^ /" sign(a;)(sin(cjt;)-^i;) 

■cos(7ra/2) / ——- du, 



■K 

which ends the proof. I 

Let a G (0, 1), then by differentiating (|15p with respect to a and taking the limit as a ^- 1~, we may 
obtain 

sign(.)|/.|log|/.HMog|.H lim ^ / signM sin(..) ^^^ ^^^^ 
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